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We calculate magnetization in magnetically doped semiconductors assuming a local exchange 
model of carrier-mediated ferromagnetic mechanism and using a number of complementary theo- 
retical approaches. In general, we find that the results of our mean-field calculations, particularly 
the dynamical mean field theory results, give excellent qualitative agreement with the experimen- 
tally observed magnetization in systems with itinerant charge carriers, such as Gai-^Mn^As with 
0.03 < x < 0.07, whereas our percolation-theory-based calculations agree well with the existing data 
in strongly insulating materials, such as Gei-^Mn^. We comment on the issue of non-mean-field 
like magnetization curves and on the observed incomplete saturation magnetization values in di- 
luted magnetic semiconductors from our theoretical perspective. In agreement with experimental 
observations, we find the carrier density to be the crucial parameter determining the magnetization 
behavior. Our calculated dependence of magnetization on external magnetic field is also in excellent 
agreement with the existing experimental data. 
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I. INTRODUCTION 

Diluted magnetic semiconductors (sometimes also re- 
ferred to as doped magnetic semiconductors - we will use 
the abbreviation "DMS" to denote both of these equiv- 
alent terminologies) have recently attracted a great deal 
of attentiorotl for their potential in combining ferromag- 
netic and semiconductor properties in a single material. 
The prototypical DMS material is Gai-^Mn^AsB (typi- 
cally x w 1 — 10%) with the Mn ions substitutionally (in 
the ideal situation) replacing Ga at the cation sites. Mn 
ions in Q&i-xM\\ x As serve a dual purpose, acting both as 
dopants (acceptors in this case) and as magnetic impu- 
rities, whose spins align at the ferromagnetic transition. 
For x s» 1 — 7%, Gai-^Mn-rAs is found to be ferromag- 
netic with the ferromagnetic transition temperature (or, 
equivalently, the Curie temperature) T c 10 — 100K. 
The optimum value of x, which corresponds to the high- 
est value of reported value of T c , is around 5%. Other 
DMS materials of current interest include Ini-JVln^ASjQ 
Gai^Mn^pJa Gei-^Mn^Ja and Gai-^Mn^Sb.Q 

In spite of a great deal of recent experimental and 
theoretical activityoo there is not yet a consensus on 
the fundamental mechanism leading to ferromagnetism in 
these systems as well as on the definitive predictive the- 
ory quantitatively describing this ferromagnetic mecha- 
nism. It is therefore important to work out detailed and 
experimentally falsifiable consequences for various pro- 
posed theoretical models and ideas. In this context, it is 
unfortunate that much of the theoretical DMS literature, 
perhaps due to the considerable technological motivation 
in creating room-temperature ferromagnetic semiconduc- 
tors for projected spintronics applications, has concen- 
trated on the calculations of T c in various DMS mate- 
rials. Such theoretical predictions of T c invariably in- 
volve tuning free parameters (e. g., the strength of the 
exchange coupling between carriers and local moments), 
whose values are often unknown a priori. This signifi- 



cantly reduces the practical importance of these predic- 
tions except perhaps in the broadest qualitative sense of 
identifying the crucial controlling parameters which de- 
termine and limit T c in DMS materials. On the other 
hand, the temperature dependence M(T) of the sponta- 
neous magnetization possesses many characteristics, such 
as concavity /convexity of the curve, value of saturation 
magnetization, critical behavior at the point of ferromag- 
netic transition etc., which cannot all be described by just 
tuning parameters of a given model. Thus study of M(T) 
has a very high potential for elucidating the physics be- 
hind DMS ferromagnetism in real systems. 

A particular issue of considerable significance in DMS 
ferromagnetism has been the non-trivial non-mean-field- 
like behavior of the spontaneous magnetization as a func- 
tion of temperature. This was. already apparent in the 
very first reported observationcl of DMS ferromagnetism 
in a IIIi-^Mn^.V material, namely Ini-^Mn^As, where 
the experimental M(T) curve exhibited an untypical out- 
wardly concave shape strikingly different from the usual 
convex M(T) behavior ^expected within the textbook 
Weiss mean-field theoryO and seen routinely in conven- 
tional ferromagnetic materials. The Ini-^Mn^As (with 
x = 0.013) system studied in Ref. ||was an insulating sys- 
tem (i. e. with resistivity increasing monotonically with 
decreasing temperature), and the insulating ferromag- 
netic DMS systems studied so far in the literature al- 
most always exhibit qualitatively similar non-mean-field- 
like concave M(T) behaviorJj 

Two of us have recently shown! that such a manifestly 
non-mean-field-like concave M(T) behavior in insulating 
DMS materials can be understood on the basis of a mag- 
netic percolation transition of bound magnetic polarons 
in the strongly, localized carrier system. Earlier numeri- 
cal simulationsEJO in the strongly localized regime had 
already indicated that the M(T) behavior of DMS fer- 
romagnets could have concave outward shapes aS|-seen 
experimentally. Very recent numerical simulations!^! in 
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the strongly localized regime have verified our polaron 
percolation picture of DMS ferromagnetism in the insu- 
lating regime. 

Even in the "metallic" (i. e. with resistivity decreas- 
ing with T below thep ferr i omami c tic traiiS |ition temper- 
ature) DMS systems^MyillMyilaiaEi such as op- 
timally doped Gao.95Mno.05As, where long range mag- 
netic ordering of the Mn magnetic moments is created 
presumably by a dilute gas of delocalized holes medi- 
ating the magnetic interaction, the experimentally ob- 
served M(T) often appears Ao be very different from 
the classic mean-field shape.B Although the "metallic" 
Ga 1 _ 2 ,Mn :I .As DMS system does not exhibit the man- 
ifestly concave temperature dependent magnetization 
seen in insulating DA^^^cmfi^.thp .observed M(T) in 
metallic DMS is ofte 1^00^ ^113130^^ almost linear 
in temperature for 0.5T C <T<T C , with a temperature- 
dependent behavior intermediate between the concave 
M(T) shapes of the localizecLtheorycl and the textbook 
convex magnetization curves.u 

We mention^,— | tkat, [ _ |[ *Ery recent annealing 
experimentMilt3taoll3Ma in G ai _ x Mn x As (with 
x = 0.05 — 0.10) demonstrate that M(T) behavior (as 
well as the value of T c ) can be strongly affected by 
annealing, and in particular, "optimal" annealing (to be 
found empirically for each sample purely by trial and 
error since the precise role of annealing in improving 
the materials quality is unknown at this stage) may 
lead to reasonably mean- field- like convex M(T) shape 
with a concomitant enhancement of T c . "Non-optimal" 
annealing, on the other hand, leads to suppression of 
T c and strongly non-mean-field-like (and non-universal) 
M(T) behavior. 

Finally, most of DMS magnetization measurements in 
the literature, particularly for the insulating DMS sys- 
tems, show a saturation magnetization (for B = 0) con- 
siderably smaller in magnitude than that expected from 
the full ordering of all the magnetic ions, indicating that 
a large fraction (sometimes as much as 90%) of the mag- 
netic ions do not contribute to the global DMS ferromag- 
netism. 

Motivated by the desire for illuminating the phys- 
ical mechanisms underlying DMS ferromagnetism, we 
theoretically consider in this paper the temperature- 
dependent magnetization M(T) in three-dimensional 
DMS systems using a number of complementary the- 
oretical approaches. The calculations we present here 
are based on the static (Weiss) molecular mean-field 
theory,Eil the dynamical mean-field theory (DMFT),E3 
and the percolation theory.c^l Of particular interest is 
the important issue of correlations between magnetic 
and transport properties of various DMS materials since 
such correlations should (and do) exist in magnetic sys- 
tems where the ferromagnetism is mediated by carri- 
ers leading to the global ordering of the dopant local 
moments. While earlier studies invariably concentrated 
on either localized-carrier ( "insulating" ) and itinerant- 
carrier ("metallic") material, we consider both these 



regimes as well as the crossover between them. We pro- 
vide detailed numerical results for our calculated magne- 
tization M(T,B), where B is the applied external mag- 
netic field and T the temperature, in various regimes 
of the system parameter space. These M(T, B) results 
should help ascertain the applicability of various theo- 
retical models to specific experimental DMS materials of 
current interest. 

This paper is organized as follows. In Sec. |l|, we de- 
scribe our model and our theoretical approaches based on 
the Weiss mean-field theory, DMFT, and the percolation 
theory, and present our numerical results for M(T, B) in 
each of the theoretical approaches providing brief discus- 
sion of our results in light of the existing experimental 
data in the literature. 
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summariz- 



We conclude in Sec. 
ing our qualitative findings and providing a critical per- 
spective on what our theoretical results on magnetization 
imply for the microscopic mechanisms underlying DMS 
ferromagnetism with the particular emphasis on the cor- 
relations between transport characteristics ( "metallic" or 
"insulating") and M(T) behavior (convex, concave, or 
"linear" ) . 



II. MODEL, THEORY, AND RESULTS 

A. The model 

We assume in this work that the fundamental mecha- 
nism underlying ferromagnetism in DMS materials (e. g. 
IIIi-zMn^V, Gei-ajMn^) is the carrier - local-moment 
(kinetic or p-d) exchange coupling, which eventually leads 
to a global ferromagnetic ordering of the impurity local 
moments {i. e. Mn) for T < T c overcoming any direct 
antiferromagnetic (super-exchange) interaction between 
the local moment spins themselves. This is certainly the 
prevalent viewpoint for DMS ferromagnetism, at least 
for Gai-^Mn^As system following the pioneering work 
of Ohno and his collaborators. It is natural to ask about 
the evidence for this belief in carrier-mediated DMS fer- 
romagnetism induced in the Mn local moments. First- 
principles band theory calculations indicate that there is 
strong p-d hybridization between the Mn d-levels and the 
valence band p-states of GaAs. This leads to strong ki- 
netic exchange coupling between hole spins and Mn spins, 
which is the basis of our model. Experimentally, there 
are a number of compelling circumstantial indications 
of Gai-^Mn^As being a carrier-mediated ferromagnetic 
material. First, there is the strong correlation (made 
even stronger by the recent annealing experiments) be- 
tween transport and magnetic properties. This, how- 
ever, is not definitively conclusive since strongly insu- 
lating Gai-^Mn^As (and Ini-^Mn^As) samples are also 
found to be ferromagnetic (albeit with much lower T c val- 
ues). Perhaps the most compelling evidence supporting 
the carrier-mediated ferromagnetic mechanism is the ob- 
served agreement (so far available only in Gai_ x Mn x As 
DMS systems) between the magnetization measured di- 
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rectly using a SQUID magnetometer and that inferred by 
analyzing the anomalous Hall effect data. Such an con- 
sistency between direct and "transport-inferred" magne- 
tization strongly suggests a carrier-mediated exchange 
mechanism underlying DMS ferromagnetism. In addi- 
tion, there are a number of experiments dealing with op- 
tical properties (again, available only for Gai^Mn^As so 
far) which indicate very strong correlation between the 
magnetic band structure and the magnetization of the 
system, leading again to the condition that the global 
ferromagnetic ordering of the Mn local moments is most 
likely induced by the hole spin polarization in DMS sys- 
tems. The related issue of whether these carriers are va- 
lence band holes or impurity band holes is more difficult 
to settle. Given the strong inherent disorder in DMS ma- 
terials and rather strong exchange coupling between hole 
and Mn spins (leading to local binding of holes and Mn 
ions), it is natural to think that all of the physics (both 
localized and delocalized) are essentially impurity band 
physics. This viewpoint, which we have adopted in our 
theories presented in this paper, has received strong sup- 
port from recent optical ex perim ents. Also our DMFT 
results, as presented in Sec. II C, clearly demonstrate the 



important role of impurity band physics in DMS mag- 
netic properties. 

Since the most intensively studied DMS materials 
of current interest are Mn-doped III-V semiconduc- 
tors IIIi-^Mn-rV (e. g. Gai_ x Mn x As, Ini-^Mn^As, 
Gai-a-Mn^P, Cai^Mn^N, Gai_ x Mn x Sb), where the 
carriers are typically holes, we will other refer to the 
semiconductor carriers as "holes" in the rest of the article 
without any loss of generality. 

So our basic model is that of a finite density n\ of mag- 
netic dopants ("impurities") interacting through a local 
exchange coupling with a finite density n c of holes in the 
host semiconductor material with n c jri\ <C 1. We as- 
sume that magnetic impurities under consideration enter 
substitutionally at the cation sites (e. g. Mn impuri- 
ties at Ga sites). Recent exp^jmrntaLa^^ studies 
of Gai_ x Mn x As have shownBBBE30li3-El that lattice 
defects may be playing an important role in determining 
magnetic and transport properties of the samples, but we 
assume in our work that these defects enter our theory 
only in determining the basic parameters of the model, 
namely, the density of magnetically active dopants n\, the 
hole density n c , and perhaps the local effective exchange 
coupling J between the holes and the magnetic impuri- 
ties, and do not include any defects into our model ex- 
plicitly. Examples of such defects which enter our model 
through its parameters rather than directly are given by 
antisite defects in the host semiconductor material (i. e. 
As at Ga sites in GaAs) and Mn interstitials (i. e. Mn 
atoms at interstitials sites rather than cation substitu- 
tional sites), which may very well be important in pro- 
viding substantial compensation in the semiconductor, 
leading to the experimental fact that the hole density n c 
is usually a small fraction of the magnetic dopant den- 
sity instead of there being a one-to-one correspondence 



between dopants and holes. The local moment density rii 
in our model is not necessarily the total Mn concentra- 
tion in the system, since the presence of Mn interstitial 
defects could lower the density of magnetically active Mn 
ions. In this work we use rij to denote local moment vol- 
ume density and x to denote the fraction of Ga atoms 
replaced by Mn dopants. 

Similarly, consistent with the spirit of our minimal 
model we also neglect all band structure effects in our 
theory, making the simplest approximation such as a sin- 
gle parabolic band with a single effective mass or a sin- 
gle simple tight-binding carrier band characterized by an 
effective band width parameter. This is not because re- 
alistic band structure effects are not of any importance 
in DMS ferromagnetism, in fact we believe that spin- 
orbit coupling in Gai^Mn^As valence band hole states 
may play a quantitative role in Gai-^Mn^As ferromag- 
netism, particularly in the relatively disorder-free "metal- 
lic" [x « 0.05) systems where the holes are likely to be 
GaAs valence band hole states with strong spin-orbit cou- 
pling. Our reasons for neglecting spin-orbit coupling and 
other band structure effects in our theory are the follow- 
ing: (1) our interest in this paper is in theoretically ex- 
ploring M(T, B) within a minimal model, which requires 
only m, n c , and J, leaving out all non-essential complica- 
tions; (2) if needed, band structure effects can be system- 
atically included in the future by appropriately extending 
the model; (3) due to the inevitable presence of strong 
exchange coupling and strong disorder in DMS materials, 
a starting point based on perfect "realistic" valence band 
hole states may be inapplicable - in fact, we believe that 
much of the DMS physics is occurring in the impurity 
band of the host semiconductor with the itinerant and 
the localized carriers being the extended and the local- 
ized hole states in the impurity band of the system and 
not the valence band states, a view point str£ 
ported by several recent experimental results.l 

The Hamiltonian of exchange interaction between 
magnetic impurities and holes we use in this paper reads 



Hm = 



J d 3 r M (Si-s(r))5(r-R, ), 



(1) 



where J, which has units of energy, is the exchange cou- 
pling between an impurity spins Sj located at Rj and 
a hole spin density s(r), and clq is the unit cell volume 
needed for proper normalization. The impurity spin S 
in Eq. ([[]) is assumed to be completely classical in the 
theory whereas the carrier spin s is treated quantum- 
mechanically. This is justified because the impurity spin 
is large, e. g. S = 5/2 for Mn in Gai-^Mn^As. 

Our mode l (ex cept for the mean-field theory considera- 
tions of Sec. |ll B|) omits the direct Mn— Mn antiferromag- 
netic exchange interaction assuming its effects to be ei- 
ther negligibly small or incorporated into the effective pa- 
rameters of the model. Actually, in the parameter range 
of interest to us (i < 1), where DMS ferromagnetism 
typically occurs, the magnetic impurities are separated 
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from each other by non-magnetic atoms, and this anti- 
ferromagnetic interaction, which rapidly decays with the 
distance, should be negligible. We also ignore any specific 
hole-hole interaction effect in our theory. These approx- 
imations are nonessential and are done in the spirit of 
identifying the minimal DMS magnetic model of inter- 
est. Both of these effects, which may be of quantitative 
importance in some situations, can be included in the 
theory by adjusting the parameters of the model or per- 
haps at the cost of introducing more unknown parameters 
characterizing these interactions. The possible effects of 
including these interactions in our calculations will be 
discussed later in the paper. 

With this introduction the full Hamiltonian is given by 



V 2 



+ 



E 



d 3 i 



]T^(r-R^t (r) ^ (r ) 



+ Ja l (S; ' ^p) S(r - R^iWV^r) 

a/3 
3 

+ /^r^^s^W^WK/a-B) , (2) 

a/3 

where the first term is the "band" Hamiltonian with m 
being the relevant effective mass, V(r) is the random po- 
tential arising from (non-magnetic) disorder, W(r — R4) 
is the Coulomb potential due to a magnetic impurity lo- 
cated at Ri, and the second term is the local exchange 
coupling between the moments of magnetic impurities 
and the carrier spins [i e. precisely the Hm defined 
through Eq. (0)], with er Q/3 being the Pauli matrices. The 
last two terms in Eq. (|^) are simply the Zeeman energies 
of the local moments and the holes respectively, g\ and 
g c are the corresponding g- factors, /is is the Bohr mag- 
neton, and B is the external magnetic field. 

The local exchange coupling can be ferromagnetic 
(J < 0) or antiferromagnetic (J > 0) without affecting 
DMS ferromagnctism [first principles band theoryEjoLill 
suggests local antiferromagnetic coupling (J > 0) for 
the holes in Gai_ x Mn x As]. Finally, we mention that 
the magnetic interaction Hamiltonian defined by Eq. (jl]) 
is sometimes—referred to as the -s^d (or s-f) exchange 
HamiltonianEj or the Zener modefij ia the literature al- 
though it was originally introducedEj by Nabarro and 
Frohlich in a slightly different context. The physics we 
are interested in is how the local exchange interaction 
defined in Eqs. (0) and (|J) could lead to global ferro- 
magnetic ordering of the impurity local moments below 
a Curie temperature T c . We mention that we choose 
S = 5/2 and s = 1/2 in all our numerical calculations 
below. 

Unfortunately none of the parameters J, n c , n\ of our 



model is directly experimentally measurable. This is why 
we have emphasized throughout this work qualitative 
behavior in temperature-dependent magnetization as a 
function of the system parameters. The carrier density 
n c is hard to measure even in metallic ferromagnetic ma- 
terials because of the problems associated with anoma- 
lous Hall effect (and the situation is obviously worse in 
strongly insulting systems). The local moment density m 
is unknown because only a fraction of the incorporated 
Mn atoms are magnetically active due to the invariable 
presence of Mn interstitials and other possible defects 
in the system. What is known about Gai-^Mn^As is 
that it is a heavily compensated system in the sense that 
the density of holes, is much less than the density of 
Mn, typically n c jn\ ~ 0.1. Much of this compensation 
most likely arises from various defects invariably present 
in low-temperature MBE grown Gai-^Mn^As. The two 
most "harmful" defects in this respect are Mn intersti- 
tials and As antisites. Both of these defects act as effec- 
tive double donors, producing two electrons each. Thus 
the holes produced by the magnetically and electrically 
active "desirable" Mn 2+ ions (sitting at substitutional 
cation sites), with each substitutional Mn ion producing 
one hole, could be heavily compensated by the defects 
leading to the existing situation in Gai-^Mn^As where 
n c jn- x <C 1. 

All aspects of Kondo physics are completely negligi- 
ble for the current problem. Kondo physics is relevant 
only in the complementary regime of high carrier den- 
sity (n c /rii S> 1), where a paramagnetic carrier ground 
state entails as the impurity spin is quenched by the free 
carriers. For itinerant delocalized carriers (i. e., metal- 
lic DMS), the indirect exchange coupling between the 
local moments (induced by carrier spin polarization) is 
precisely the RKKY interaction. The relevance (or per- 
haps even, the dominance) of RKKY physics (leading 
to a magnetic ground state) over the Kondo physics in 
the low carrier density limit of the Kondo lattice sys- 
tem has occasionally been mentioned in the Kondo-effect 
literatures 

The crucial element of RKKY physics playing a key 
role in DMS ferromagnetism is the relatively low values of 
carrier density in these materials leading to kp^ni <C 1 
(where kp is the Fermi wavevector associated with the 
carrier density n c and ^/n{ is the characteristic inter- 
impurity separation) so that the RKKY interaction is 
essentially always ferromagnetic, and the RKKY spin- 
glass type behavior predominant in disordered magnetic 
metallic systems may not arise here since the interaction 
is mostly ferromagnetic avoiding effects of frustration. 

In principle one could start from the general Hamilto- 
nian (^), and try to develop a theory for carrier localiza- 
tion and magnetism on an equal footing. Such an ambi- 
tious attempt would be essentially futile, due to the enor- 
mous complexity of the problem, which would require 
both disorder and exchange interaction to be treated non- 
perturbatively. We therefore adopt the reasonable empir- 
ical approach of building into our basic model the metal- 
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lie (itinerant carriers) or the insulating (localized carri- 
ers) nature of the system, and develop separate comple- 
mentary theoretical approaches for the two situations in 
order to compare and contrast the nature of DMS ferro- 
magnetism in metallic and insulating systems, both with 
the local exchange magnetic Hamiltonian (fy) coupling 
the impurity moments to carrier spins. 

We use complementary theoretical approaches 
(degenerate- and non-degenerate-carrier Weiss static 
mean- field theory, dynamical mean field theory (DMFT), 
and percolation theory) in this paper; two of which 
apply to the limiting cases of extended "metallic" system 
(degenerate-carrier mean-field theory and DMFT) and 
the other two to the strongly localized "insulating" 
system (non-degenerate-carrier mean-field theory and 
percolation theory). In principle, DMFT can interpolate 
smoothly between the systems with extended and 
localized carriers, but we have neglected localization 
effects in our DMFT calculations carried out so far. Our 
results for temperature- and magnetic-field-dependent 
magnetization M(T, B) exhibit qualitative behavior 
very similar to that seen in experiments which, given 
the minimal nature of our theoretical model, is all we 
can expect of the theory. In our view, at this early 
stage of the development of DMS materials and their 
physical understanding (and given the metastable and 
fragile complex nature of the DMS systems), where all 
the materials details (e. g. defects and disorder of many 
possible types) qualitatively affecting the system magne- 
tization may not even be known at the present time, it is 
premature to demand (or impose, for example, by tuning 
the parameters of the theoretical model) quantitative 
agreement between theory and experiment. In this re- 
spect we strongly disagree with the statements asserting 
that DMS ferromagnetism is a well-understood problem 
based on a model of free valence band holes interacting 
with Mn local moments. The excellent qualitative 
agreement between our calculated magnetization results 
in different theoretical approaches and the experimental 
magnetization data in the existing literature provides 
useful insight into the possible magnetic mechanisms 
underlying "metallic" and "insulating" DMS systems. 



B. The Static Mean-Field Theory 

The basic idea underlying the static mean-field the- 
ory, as applied to ferromagnetism in DMS, is to repre- 
sent action of all impurity/hole spins upon a given impu- 
rity/hole spin as an effective "mean field," whose value is 
determined by the average values of the spins acting upon 
this given spin. The resulting equations for the spins of 
impurities and holes are to be solved self-consistently, fi- 
nally yielding the equilibrium magnetization at a given 
temperature. 

The difference between the mean-field theory consid- 
ered in this Section and the canonical Weiss mean-field 
model arises from the existence of two interacting species 



of spins, those of holes and of impurities in a DMS sys- 
tem. As a result, we have two effective fields - one deter- 
mined by the average value of a hole spin and the other 
determined by the average value of an impurity spin. In 
the framework of Hamiltonian (Q), the effective field act- 
ing upon holes has contributions coming from magnetic 
impurities and from the external magnetic field B, 



B 



(c) _ 



1 



(JalmiSz)) + B 



(3) 



where the direction of the z axis is chosen to coincide with 
the direction of applied magnetic field B, or, in the case of 
-6 = 0, with the direction of spontaneous magnetization 
of impurities. The effective field acting upon impurities is 
a sum of contributions from holes and the external mag- 
netic field. Relative simplicity of the static mean-field 
theory allows us to account also for possible direct an- 
tiferromagnetic interaction between magnetic impurities 
by adding the following term to the Hamiltonian 



jk 



(4) 



which yields one more contribution to the effective field 
acting upon magnetic impurities: 

BS = — {z AF J AF (S z ) + Ja 3 n c ( Sz )) +B, (5) 

where z AF is the effective number of surrounding impu- 
rities a given impurity interacts with. 

The response of the impurity spin to this effective field 
B^ff is given by 



(S z ) = SB S 



k B T 



(6) 



where 



B s {x) 



2s 



2s 



1 2s 
■ coth ■ 



2s 



— coth —x (7) 
2s 2s w 



is the Brillouin function. The magnetic response of 

(c) 

hole spins to effective field B^ produced by impuri- 
ties strongly depends on whether the hole gas is de- 
generate or not. The two complementary cases of non- 
degenerate a nd de gene rate h oles will be considered be- 
low in Sees. II B I and II B 2 respectively. We mention 
that the non-degenerate (degenerate) situation applies 
primarily to the insulating ("metallic") DMS systems. 



1. Non-degenerate holes 

The case of non-degenerate holes, when the hole spin 
distribution is not affected by the Pauli exclusion princi- 
ple, corresponds to two physical situations. The first one 
occurs when the holes are localized with strong on-site 
repulsion, so there is only one hole at each localization 
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center. This scenario is relevant to DMS with l ocali zed 
carriers, which will also be considered in Sec. IIP us- 
ing the percolation theory. The second situation takes 
place in a gas of delocalized holes when the temperature 
is higher than the Fermi energy. 

When the Pauli exclusion principle plays no role in the 
spin distribution of electrons, the latter is determined 
by Boltzmann statistics, and the average hole spin, as 

(c) 

determined by the effective mean-field £r ff is given by 



SB S 



k R T 



(8) 



similarly to Eq. (p), with g c being the hole ^-factor. 



Combining Eqs. (§), (§), (§), and (§, we obtain the 
self-consistent equation for (S z ): 



S 



= Bs 




t V m 'V (s + i)(5 + l) 

6STi (S z ) giUBB 




(S + 1)T S k B T 



(9) 



where 



k B T c0 = -Jaly/nJi[y/S(S + l)s(s + 1) 



AF 7 AF 



fcsTi = -S(S + l)z^J 



(10) 



(11) 



and J is assumed to be negative (antiferromagnetic in- 
teraction between impurities and holes). Similarly to the 
text-book mean-field theory, Eq. (0) has a non-trivial so- 
lution in the absence of the external magnetic field only 
if temperature T is below a certain value, which is the 
ferromagnetic transition temperature T c . Using the ex- 
pansion for the Brillouin function 



Bs(x)\ 3 



1 



we arrive at the following expression 



x + 0(x 3 ) , 



T« = */23> + Z?-Ti. 



(12) 



(13) 



In the absence of antiferromagnetic interaction be- 
tween the impurities, T± = so the ferromagnetic transi- 
tion temperature equals T c0 . Equation (ph in the absence 
of the external magnetic field and at J = reduces to 




the average hole spin is still given by Eqs. (||) and (||) 
with B = 0. The solution of Eq. (14) can be found 



numerically. The resulting magnetization curves for im- 
purities and holes for several values of n c /nj are shown 
in Fig. |l|. Note that specific values of J, ni, and n c are 
not of any relevance here - only the ratio n c jri\ is the im- 
portant tuning parameter in determining magnetization 
as a function of T/T c . The most important salient fea- 
ture of Fig. is the highly "non-mean-field-like" concave 
magnetization behavior for low values (< 0.2) of n c jn\. 
The reason for such behavior is that if we have one hole 
per many impurities, the effective field £T ff acting on 

holes is much stronger than its counterpart acting 
on impurities. As a result, the hole magnetization grows 
as 



iff: 
•s 



/ T (n)_ T 



T 



for T < T c o and reaches unity at some temperature of 
the order of T c o, while the impurity magnetization is 
still much less than unity. As the temperature is get- 
ting lower, the impurity magnetization grows as 




(,s + l)(5+l) 



and approaches unity only at 



Finite antiferromagnetic coupling J AF suppresses T c , 
as one can easily see from Eq. (|30|). In Fig. || we show 
the dopant magnetization, as given by the (numerical) 
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FIG. 1: (a) Dopant magnetization M/M = (S g )/S and (b) 
hole magnetization m/mo = (s z )/s with S = 5/2 and s = 1/2 
from the non-degenerate-hole model for various density ratio 
{n c /m = 1.0, 0.5, 0.2, 0.1, 0.05). 



solution of Eq. (0) for B = 0, for various values of Ti/T c q 
for two typical values of n c jri\ — 0.5 and 0.1. 

In discussing the results shown in Fig. || we first note 
that the actual direct Mn-Mn antiferromagnetic coupling 
J AF is expected to be very small in Gai-^Mn^As for 
low values of x, i.e. for relatively large Mn-Mn sepa- 
ration, since the direct antiferromagnetic coupling falls 
off exponentially with inter-atomic distance. For larger 
values of x, however, effects of antiferromagnetic Mn- 
Mn coupling may very well be important in determining 
the DMS magnetic behavior. In general, finite J AF sup- 
presses both T c and M(T) as one would expect. In par- 
ticular, the zero-temperature magnetization M(T — » 0) 
could be strongly suppressed far below the saturation 
magnetization by the finite AF coupling, particularly for 
lower hole densities, 



s 



T^0 




An interesting feature of Fig. || is that at low carrier 
densities larger values of J AF may actually restore (albeit 
with strongly suppressed value of saturation magnetiza- 
tion) the convex Af(T) shape (e. g. the z AF |J AF |/J = 
5.0 x 10~ 3 curve for n c jni = 0.1 in Fig. ||(b)). For higher 
hole densities, however, M(T) becomes more concave for 
J AF £ 0. 

The magnetic susceptibility of the system is essentially 



FIG. 2: Dopant magnetization with both J AF and J from 
the non-degenerate- hole model with S = 5/2 and s = 1/2 for 
various coupling constant ratio (z AF J AF /J = 0.0, 0.5, 1.0, 
2.5, 5.0xl0 -3 , from the top), (a) for density ratio n c /rii = 0.5 
and (b) n c /rii — 0.1. 



that of impurities (since n c /rii -C 1 and S > s), and 
above the critical temperature is given by 



X(T) = 



_ d{giH B {S z )) 



Xo 



dB 
T 2 Jal 




(16) 



where T* = T x + ^JT'( + T c 2 , and X o = ni{g i( i B ) 2 S{S + 
1)/3T is the paramagnetic susceptibility of bare impuri- 
ties. 

Finally, we show in Fig. || the effect of an external mag- 
netic field by showing M(T, B) for fixed parameter values 
J, S, s, and for various density ratio n c jn\. In the main 
figure we show M(T, B) as a function of temperature for 
a magnetic field value, B = 2.6T (solid lines), and in the 
inset we show M(T, B) as a function of the external mag- 
netic field for different temperatures. In inset we show 
M(B,T C ) oc i? 1 / 3 as expected in the mean-field theory. 
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T=0.95, 1.0, 1.05 T c 




FIG. 3: External magnetic field and temperature dependence 
of dopant magnetization for the non-degenerate-hole model 
for various density ratio (n c /n; = 1.0, 0.5, 0.2, 0.1, 0.05, from 
the top). Solid (dashed) lines indicate the results for a fixed 
external magnetic field B = 2.6T (0T). The parameters of 
Fig. 1 are used. In the inset the magnetization curves as 
a function of magnetic field are given for fixed temperatures 
(T =0.95, 1.0, 1.05T C , from the top). Solid (dashed) lines 
represent the results for n c /ni = 1.0 (n c /ni = 0.1). 



2. Degenerate hole gas 

For metallic DMS systems where ferromagnetic T c is 
optimum (at least for Gai_ x Mn x As), the carrier system 
is typically delocalized, with the Fermi energy substan- 
tially exceeding the ferromagnetic transiti on te mpera- 
ture. In such systems, the formalism of Sec. II B 1 , devel- 
oped for the non-degenerate hole system is not applicable 
anymore, and we have to use the generic expression 



s 



1 1 

n c 2 



def(e) D{e + g cf i B sB^) 

D(e - g B HB8B®j\ (17) 



for the hole polarization, where D{e) is the hole density 
of states. 

(c) 

If the effective field B^g acting on carriers is weak, we 
can expand the density of states up to the first order in 
the effective field to obtain 



S 



-g cf iBsB^ f def(e)D'(e) 



(18) 



For very low temperatures, ksT <C Ep (in general, this 
condition is satisfied in metallic Gai-^Mn^As systems 
where typically ksTc/Ep < 0.1 noting that both T c and 
Ep decrease with lowering the hole density n c ), we get 

— = -9^bsB^D'(E f ) . (19) 



Linearizing Eq. (||) and using Eqs. (0) and (||), we arrive 
at 



7 1 i d) = - 2Ti , 



where 



, (d) S(S + 1) 



L c0 



^{Jal) 2 D{E F ) ni 



(20) 



(21) 



and T\ is given by Eq. (|1 1|) . 

In the absence of antiferromagnetic interaction be- 
tween magnetic impurities, J AF = 0, the ferromagnetic 
transition temperature equals T c q. As one can see from 
Eq. (|2l|), T c0 cx J 2 nl' 3 in the Weiss mean-field the- 

1/3 

ory for degenerate holes (since D(Ep) oc n c for a 
three-dimensional degenerate electron gas) , in con trast to 
the non-degenerate case considered in Sec. |ll B 1 , where 
T c q oc Jr^J 2 . 

When the hole density is very small, the hole Fermi en- 
ergy Ep may be comparable to the effective magnetic en- 
ergy g c (J-BsB^Q . In this case we can not expand Eq. ( p"7j ) 

(c) 

with respect to the effective field B^ s ', and the hole 
magnetization must be obtained by directly integrating 
Eq. (g§). 

Before presenting our numerical results for the cal- 
culated temperature-dependent DMS magnetization for 
the degenerate-carrier mean-field theory, we mention 
that equation ( pp| ) for the ferromagnetic transition 
temperature in the Zener model was first derived in 
Rcf . |35| morc_than forty years ago and has recently been 
rcdiscovcrcaOEZI in the context of DMS ferromagnetism. 

In Fig. |] we show our calculated impurity and car- 
rier magnetization for the degenerate-carrier Weiss mean- 
field theory using the same parameters (typical for 
Gai^aMn^As) as for the corresponding non-degenerate 
case shown in Fig. pi. In general, except for the lowest 
hole density with n\Jn c = 0.05, the impurity magnetiza- 
tion is convex, but looks quite different from the classic 
text book Weiss form except perhaps at very high carrier 
densities when rii/n c « 1. In particular, for the realis- 
tic value of n c jn\ = 0.1 (which is thought to apply to 
many Gai-^Mn^As samples), M(T) in Fig. |] has the 
non-mean-field-like straight line shape at lower temper- 
atures. For higher values of hole density, n Q jn\ =0.2 
and 0.5, the magnetization is much closer to the classical 
Brillouin shape whereas for lower carrier densities (e. g., 
n c /rii — 0.05) even this metallic DMS system starts ex- 
hibiting the concave magnetization curve typical of the 
i nsulat ing DMS systems discussed in the last sub-section 
( jHBl| ) of this article. 

The origin of the concave shape of the magnetization 
curve is exactly th e sam e as for the non-degenerate case 
discussed in Sec. II Bl, namely, the magnetization of 
holes saturates at T < T c q, much earlier than magne- 
tization of impurities, which, upon saturation of the hole 
magnetization grows as given by Eq. (|l4|), which repre- 
sents a concave curve provided n c /n; <C 1. Thus, the 
concave magnetization behavior may be generic to the 
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low carrier density limit of DMS ferromagnets, indepen- 
dent of whether they are metallic or insulating although 
the highly concave magnetization curves of Fig. |](a) are 
clearly much more typical of insulating DMS systems 
than the metallic ones. The carrier magnetization re- 
sults presented in Fig. [|(b) are very similar to textbook 
convex Weiss magnetization behavior. We note that 
the degenerate-carrier mean-field theory results shown 
in Fig. ^a) are qualitatively very similar to the exper- 
imentally measured temperature dependent magnetiza- 
tion data in metallic Gai-^Mn^As systems. In particu- 
lar, recent annealing experiments, where annealing leads 
to better metallicity (e. g., higher conductivity) by virtue 
of increasing the hole density, show qualitative trends 
strikingly similar to the results of Fig. |(a). 

Next we include the direct antiferromagnetic coupling 
(||) between the local moments in our consideration. 
With J AF ^ 0, the critical temperature is given by 
Eq. (|2^). The calculated impurity magnetization in the 
presence of the antiferromagnetic coupling is shown in 
Fig. [|. Comparing with the non-degenerate-hole model 
we find that the degenerate hole case is more quanti- 
tatively sensitive to the antiferromagnetic coupling al- 
though the qualitative effect of a finite J AF in both metal- 
lic and insulating DMS materials is basically the same, 
namely suppression of the transition temperature and the 
magnetization. 

Similar to the results shown for the insulating systems, 
we present in Fig. [5] calculated impurity magnetization 
for two values of the ratio n c /n; = 0.5 [Fig. ||(a)] and 0.1 
[Fig. |(b)]. 

In Fig. |6| the calculated external magnetic field and 
temperature dependence of the dopant magnetization is 
shown for fixed parameter values J, S, s, and for var- 
ious density ratio n c /rii. In the main figure we show 
M(T, B) as a function of temperature for a magnetic 
field value, B = 2.6T (solid lines), and in the inset 
we show M(T, B) as a function of the external field for 
fixed temperatures. Our calculated M(T, B) behavior 
is roughly qualitatively similar to experimental observa- 
tions in G&i- x yhi x As systems in the "metallic" regime. 

The magnetic susceptibility of the syste m is essentially 
that of impurities, similarly to Sec. II Bl, Eq. (|l6|), and 
above the critical temperature is given by 



T, 



(d) 



X(T) = 



X0 + rM 9i9cf * B 



T, 



(d) 



(22) 



T 



where xo = ^(giUs) S(S + 1)/3T is the paramagnetic 
susceptibility of bare impurities. 

The Weiss molecular mean- field theory results (for lo- 
calized car riers in Sec. [I B 1 and for delocalized carri- 
ers m Sec. |IIB2| ) presented above for DMS magnetiza- 
tion qualitatively agree very well with the existing ex- 
perimental data. In particular, the basic trend of our 
results shown in Figs. that the spontaneous magne- 
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values: S = 5/2, s = 1/2, n\ = 
x — 0.05, and coupling constant J 



FIG. 4: (a) Dopant magnetizations from the degenerate-hole 
model for various value of the density ratio (n c /ni = 0.5, 0.2, 
0.1, 0.05, from the top). Here we use the fixed parameter 

10 21 cm -3 , m = 0.5m e , 
= 3.0 eV. (b) Hole magne- 
tization with the same parameters as in (a). Note m(T — 0) 
is the magnetization at T = 0, which may not be equal to 
m = g c fJ.Bsn c . 



tization is strongly suppressed (perhaps even into a very 
unusual outwardly concave shape) at low carrier densities 
and for more insulating systems whereas at higher carrier 
densities and for more metallic systems the magnetiza- 
tion has the usual convex textbook shape, is in excel- 
lent agreement with experiment. Our degenerate-carrier 
mean- field results (see Fig. 3(a)) also reproduce the al- 
most linear magnetization curves seen in GaMnAs for 
intermediate carrier densities. But obviously one needs 
to go beyond the Weiss mean-field theory for a deeper 
(and more quantitative) understanding of DMS magne- 
tization, if for no other reason than to validate (or to as- 
certain the regime of validity of) the simple static mean- 
field theory. At low carrier density the mean-field theory 
should work better since on the average there are many 
Mn local moments between any two holes, but typically 
the number of Mn atoms per hole is around 3-4 so the 
quantitative applicability of MFT is questionable. Note 
that the hole-hole interaction neglected in our calcula- 
tions may be included in a crude approximate fashion 
by incorporating a Stoner enhancement of the carrier 
susceptibility, though in general the strength of this en- 
hancement is unknown. 

In the next two sections we go beyond the Weiss static 
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FIG. 5: Dopant magnetization with both J 7^ from the 
degenerate-hole model for various values of the coupling con- 
stant ratio (a) z AF J AF /J = 0.0, 0.5, 1.0, 1.5, 2.5 xlO" 3 (from 
the top) and for density ratio n c /ni = 0.5, (b) z AF J AF /J = 
0.0, 0.5, 1.0, 1.5, 1.75X10" 3 (from the top) and n c /m = 0.1. 
The parameter values of Fig. ^ are used. 



mean-field theory and develop two more sophisticated ap- 
proximation schemes to theoreticall y stu dy DMS magne- 
tization. These arc D MFT (section II C ) and the perco- 
lation theory (section |ll D|) . 



C. The Dynamical Mean-Field Theory (DMFT) 

In this section we use a recently developed non- 
perturbatuje method, the "dynamical mean-field theory" 
(DMFT),l3 to calculate magnetization for the minimal 
model (Eq. (||)) of dilute magnetic semiconductors. The 
DMFT has been recently applied to the DMS system to 
calculate the magnetic transition temperature and the 
optical conduct ivityEE DMFT is essentially a lattice 
quantum version of the Weiss mean field theory where the 
appropriate density of states (including impurity band 
formation) along with temporal fluctuations are incorpo- 
rated within an effective local field theory. 

We model the Gai-^Mn^As system as a lattice of sites, 
which are randomly nonmagnetic (with probability 1 — x) 
or magnetic (with probability x), where x now indicates 
the relative concentration (i.e. per Ga site) of active Mn 
local moments in GaMnAs. The DMFT approximation 
amounts to assuming that the self energy is local or mo- 
mentum independent, E(p, zcj„) — > Yl(iu) n ), and then all 
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FIG. 6: External magnetic field and temperature dependence 
of dopant magnetization from the degenerate-hole model for 
various values of the density ratio (n c /m = 0.5, 0.2, 0.1, 0.05, 
from the top). Solid (dashed) lines indicate the results for a 
fixed external magnetic field B — 2.6T (0T). The parameters 
in Fig. ^| are used. In the inset the magnetization curves as 
a function of magnetic field are given for fixed temperatures 
(T =0.95, 1.0, 1.05T c (d) , from the top). Solid (dashed) lines 
represent the results for n c /ui = 0.5 (n c /rii = 0.05). 



of the relevant physics may be determined from the local 
(momentum-integrated) Green function defined by 



Gloc (ioJn 



4 



d 3 p 



1 



(27r) 3 iui n + [i + ha - e(p) - S CT (iw„) ' 

(23) 

where we have normalized the momentum integral to the 
volume of the unit cell a§, and /1 is the chemical potential 
and h the external magnetic field. G\ oc is in general a 
matrix in spin and band indices and depends on whether 
one is considering a magnetic (a) or non-magnetic (b) site. 
Since Gi oc is a local function, it is the solution of a local 
problem specified by a mean- field function which is 
related to the partition function Z\ oc = J dS cxp(— 5i oc ) 
with action 



S 



f (S • cr a p)ct a (x)c u0 (x) , (24) 

ua.fi 



on the a (magnetic) site and 

S^ C =gl af AT-T')ci a {T)c a p(T>), (25) 

on the non-magnetic (b) site. The local Green function 
is calculated exactly as 



G 



k>c(iw n ) = ^(c 



JS • 



(26) 
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where the thermal average (...) is taken with respect to 
the orientation of the local spin S. The a-site mean- field 
function can be written as g^ a n = a + aim • er Q/ g with 
m the magnetization direction and ai vanishing in the 
paramagnetic state. It is specified by the condition that 
the local Green function computed from Zi oc , namely 
Shi Z\ oc /SgQ = (<7g — E) - is identical to the local Green 
function computed by performing the momentum inte- 
gral using the same self energy. 

The form of the dispersion given in full Hamiltonian 
Eq. (g) applies only near the band edges. It is nec- 
essary for the method to impose a momentum cutoff, 
arising physically from the carrier band- width. We im- 
pose the cutoff by assuming a semicircular density of 
states D(e) =a%f -^S(e - e pa ) = V4i 2 - £ 2 /2nt with 

t = (2-7T) 2 / 3 /ma,Q. The parameter t is chosen to correctly 
reproduce the band edge density of states. Other choices 
of upper cutoff would lead to numerically similar results. 
This choice of cutoff corresponds to a Bethe lattice in in- 
finite dimensions. Other (perhaps more realistic) choices 
for the density of states would give magnetization results 
qualitatively similar to our results since the band edge 
density of states has the correct physical behavior in our 
model. For this N(e) the self consistent equation for g 
obeys the equation 

g g(u;) = g£(w) = u + »-xt 2 (( s %(cj) + js-tr afi y 

(l-xtfi&u)- 1 ), (27) 

where the angular brackets denote averages performed in 
the ensemble defined by the appropriate Z\ oc . 

Within this approximation the normalized magnetiza- 
tion M(T) of the local moments is given by 



M{T) = / (dS) ■ S 



,ex p(-,Sj oc ) 

Z\oc 



(28) 



As the temperature is increased, the spins disorder 
and eventually the magnetic transition temperature is 
reached. Above this temperature, go is spin- independent. 
By linearizing the equation in the magnetic part of go 
with respect to a± we may obtain the ferromagnetic tran- 
sition temperature T c . The details on the calculation of 
T c are given in Ref. |2^. The critical temperature T c and 
magnetization of the local moment depend crucially on 
J/t, x, and carrier density n c . Note that in the DMFT 
calculation it is more convenient to express our results 
in terms of n, the relative concentration of active local 
moments rather than m, the absolute impurity density. 

Before we show our calculated magnetization we de- 
scribe the density of state (DOS) of dynamical mean- 
field calculations applying to simple semi-circle models. 
In Fig. fj] we show the DMFT density of states for major- 
ity spin near the band edge corresponding to the disor- 
dered spin state T > T C) and the inset shows the T = 
ferromagnetic state. The evolution of the energy (ui) de- 
pendent DOS is shown as the carrier-spin coupling J is 
increased from zero; note that the method works equally 
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FIG. 7: The calculated density of states for dynamical mean 
field calculations applied to the semicircular model. Shown 
is the evolution of majority spin DOS for various carrier-spin 
coupling J in the disordered spin T > T c state. Inset shows 
the DOS for the ordered ferromagnetic state at T = 0. For 
the large coupling constant J > J c we find a well separated 
impurity band below the main band. 



well in the ordered T = state and the disordered spin 
T > T c case, and predicts the formation of a spin po- 
larized impurity band for J > t. For small J we see 
the expected band shift proportional to x J. For J > J c 
an impurity band centered at ~ — J and containing x 
states is seen to split off from the main band, where the 
critical coupling J c ~ t. In the DMFT calculations we 
parametrized the exchange coupling J in terms of the 
band width parameter t, subsuming the unit cell volume 
Og implicitly. 

In Fig. U we show the normalized magnetization of the 
local moments as a function of temperature for differ- 
ent values of J =1.0, 1.5, 2.0t and x = 0.05, and for 
various hole densities, n c /rii. For the small coupling con- 
stant J — t the impurity band is not formed, but for 
J =1.5, 2.0£ we have a spin polarized impurity band. 
For relatively high density (n c /rii — 0.4, or n = 0.02) 
the calculated magnetization looks similar to the Weiss 
mean-field results. But for low density (n c /ni = 0.04) 
we have a linear M(T) in the intermediate temperature 
range. Near the critical temperature T c the critical be- 
havior of the magnetization for all density is given by 
M(T) cx (T c — T) 1 / 2 . For different exchange couplings 
we have similar results (i e., linear behavior at low den- 
sities and intermediate temperature ranges). 

In discussing our DMFT magnetization calculations we 
note that the DMFT results shown in Fig. || are qualita- 
tively roughly similar to the delocalized static (degener- 
ate) mean-field-theory results we obtained in Sec. |[I B 1 . 
For example, the results of Fig. || approximately resem- 
ble those shown in Fig. 0(a) except that the DMFT 



12 




0.0 0.2 0.4 0.6 

T/T e 

l.o ' 1 ~" ~ 



8 1.0 



(b) 




0.0 0.2 0.4 0.6 0.8 

T/T. 

1.0 



1.0 




0.0 0.2 0.4 0.6 0.8 1.0 

T/T c 

FIG. 8: The normalized DMFT impurity magnetization as a 
function of temperature for (a) J = 1.5t, (b) J = l.Ot, and (c) 
J = 2.0t, for x = 0.05 and for n c /m =0.4, 0.2, 0.1, 0.04 (from 
the top). The dashed line in (a) represents the magnetization 
calculated for the simple Weiss mean field theory for the local 
moment spin S — 5/2. 



results of Fig. || are in much better agreement with the 
magnetization measurements in metallic GaMnAs sys- 
tems in the sense that the linear behavior of M(T) for 
lower T (with almost a kink just below T c as can be 
seen in Fig. ||) is much more pronounced (as in the ex- 
perimental data) than our delocalized mean-field-theory 




T/T c 



FIG. 9: The external magnetic field effects on the magneti- 
zation as a function of temperature for fixed parameter values 
J = 1.5t, x — 0.05, and n c /rii = 0.2. The curves correspond 
to B=0.002, 0.001, 0.0005, 0.0001, O.Ot (from the top). In- 
set shows the magnetization as a function of external field at 
T = T c . 



results shown in Fig. ^(a). This is both gratifying and 
expected because DMFT is a substantial improvement 
on the Weiss MFT as it incorporates the physics of spin- 
split impurity bands through the appropriate quantum 
self-ener gy c orrections not included in the simple MFT 
of Sec. II B. In fact, for very low carrier densities we 



obtain outwardly concave M(T) in our DMFT calcula- 
tions (as we expect to do in the strongly nondegenerate 
limit), but the computational convergence in our DMFT 
numerical calculations is rather poor in this regime of 
unrealistically low (n c /rii < 0.01) carrier densities, and 
therefore we refrain from showing these results. One can, 
however, detect very slight concavity in the lowest car- 
rier density (n c /rii = 0.04) M(T) results shown in Fig. 
|8|. We point out that the critical magnetic properties in 
DMFT are the same as in the static MFT, and therefore 
all the DMFT critical exponents are equal to those in the 
Weiss MFT. 

In particular, it should be possible to obtain the M (T ) 

H) 



IIBl 



behavior for the localized carrier case (c/. Sec. 
also from DMFT by incorporating impurity band local 
ization in the DMFT formalism. Our current theory does 
not include localization, and the impurity band (or va- 
lence band) carriers in our DMFT calculations are all 
delocalized metallic carriers. First principles band the- 
ory calculations indicate that the actual exchange cou- 
pling in GaMnAs may be close to critical J c , and as such 
impurity band physics may be quite important for un- 
derstanding DMS magnetization. To incorporate physics 
of localization one needs to include disorder effects (in- 
variably present in real DMS systems) in th e mo del. All 
our MFT calculations (both DMFT in Sec. |lC| and the 
static MFT of Sec. HB) are done in the virtual crystal 



approximation where effective field is averaged appropri- 
ately leaving out random disorder effects explicitly. In 
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the next section we explicitly incorporate disorder in the 
theory by developing a percolation theory approach to 
DMS magnetization for the strongly localized insulating 
systems. 

The calculated magnetization as a function of temper- 
ature is shown in Fig. H for various external magnetic field 
values. The inset shows the magnetization as a function 
of external field at T = T c . Our calculated M(T, B) be- 
havior is roughly qualitatively similar to the mean-field 
results in section [I B 2 . At the critical temperature T c 
M (T c , B) shows a mean-field behavior, M (T c) B) cx B 1 / 3 . 
The DMFT results shown in Fig. ^ are qualitatively sim- 
ilar to DMS experimental results. 



a 

o 
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D. Magnetization in the percolation formalism 



Our recently developed percolation theoryO applies 
strictly in the regime of strongly localized holes where 
the dynamical mean-field theory for delocalized carriers 
described in Sec. [I C has little validity. These two theo- 



ries, mean-field theory and percolation theory, are there- 
fore complementary. Interestingly, however, the pertur- 
bation theory and the dynamical mean-field theory are 
not mutually exclusive in spite of their regimes of validity 
being different, and in particular, a significant aspect of 
our percolation approach is its ability to reproduce qual- 
itatively the mean-field theory results of the last section 
both for T c and M (T) . 

The percolation theory assumes the same carrier- 
mcdiatcd fcrromagnetism model of Sec. II B, but now 
the carriers are pinned down with the localization radius 
05. The localized carriers are therefore taken to be in 
the impurity band and disorder, completely neglected in 
the mean- field theory of Sec. II B and [IC, now plays 
a key role in the carrier localization. The mean-field 
theory and the perturbation theory are therefore com- 
plementary in the sense that one (the mean-field theory) 
completely neglects disorder, and the other (the percola- 
tion theory) includes disorder at a very fundamental level 
i. e. by starting from the picture of localized carriers in 
a strongly disordered system. The magnetic impurities 
in this case are assumed to be completely randomly dis- 
tributed in the host semiconductor lattice in contrast to 
the mean-field case where the carrier states are free and 
the disorder is neglected. 

As we have demonstrated in Ref. the problem of fer- 
romagnetic transition in a system of bound magnetic po- 
larons can be rigorously reduced to the problem otpver- 
lapping spheres studied in the percolation theory.E3 The 
latter problem studies spheres of the same radius r ran- 
domly placed in space (three-dimensional in our case) 
with some concentration n. Overlapping spheres form 
"clusters;" as the sphere radius r becomes larger, more 
and more spheres join into clusters, the clusters coalesce, 
and finally at some critical value of the sphere radius an 
infinite cluster spanning the whole sample appears. This 
problem has only one dimensionless parameter, r 3 n, and 



FIG. 10: Curie temperature T c as a function of the dimension- 
less parameter a%n c . At a%n c < 1, T c is given by Eq. (po[). At 
a%n c > 1, Eq. (BOj) (being beyond limits of its applicability ) 
predicts decline of T c (dotted line); in reality, T c grows mono- 
tonically with a%n c (solid line), though its exact behavior is 
unknown. 



therefore can be easily studied by means of Monte-Carlo 
simulations. 

Each sphere of the overlapping spheres problem corre- 
sponds to a bound magnetic polaron, which is a complex 
formed by one localized hole and many magnetic impu- 
rities with their spins polarized by the exchange interac- 
tion with the hole spin. The concentration n of spheres is 
therefore equal to the concentration n c of localized holes. 
The expression for the effective polaron radius is not triv- 
ial and has been found in our earlier work.EI The result- 
ing formal relation between the physical parameters of 
the system under consideration and the Q-ply parameter 
of the overlapping spheres problem reads u 



r 3 n 



0.86+ (a 3 R n c 



T 



(29) 



Here 0.64 w 0.86 3 is the critical value of the parame- 
ter r 3 n at which the infinite cluster appears, and T c is 
the Curie temperature of the ferromagnetic system under 
consideration, derived in Ref. O, 



T c ~ sSJ 




exp 



0.86 



The limit of applicability of Eq. 



(30) 

( pOj ) is determined by 
the condition a 3 B n c <C 1. The dependence of the Curie 
temperature on the hole concentration in a wider domain 
of values of parameter a 3 B n c is shown schematically in 
Fig. [It]. 

Since the magnetic characteristics of the sample are 
mostly due to magnetic impurities rather than holes, the 
quantity of interest is the number of magnetic impurities 
in a cluster, which is proportional to the cluster volume. 
The magnetic properties of the system can be expressed 
in terms of the following quantities, which can be easily 
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T/T n 
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FIG. 11: Fractions of volume taken by the infinite cluster 
(solid line) and finite clusters (dashed line) and the fraction 
of volume that do not belong to any cluster or sphere (dotted 
line) for a 3 B n c = 10 -3 (Monte-Carlo simulation) 



found: (i) concentration V(vn;r 3 n)dv of clusters with 
volume between v and v + dv, (ii) the fraction of vol- 
ume Poo(r 3 n) taken by the infinite cluster, and (Hi) the 
fraction Po(r 3 n) which does not belong to any sphere 
or cluster of spheres. Clearly, these quantities obey the 
relation 

/>oo 

Poc{rn^)+ / V(vn;rn^) v dv + P (rni) = 1 . 
Jo 

Figure |ll| shows the behavior of the three terms of the 
above equation as a function of temperature T, with r 3 n 

= 1Q- 3 . 



related to T/T c by Eq. (|2Sj) and a 3 B r 

Having found these quantities, one can easily find the 
magnetic properties of the system. 

The spins of non-connected clusters are not correlated 
and average out when we calculate the spontaneous mag- 
netization of the whole sample. The only non-vanishing 
contribution comes from the infinite cluster. The total 
magnetic moment of the sample per unit volume is there- 
fore 



T 

M =n i SP 00 [ — ,a%n ( 



(31) 



Here and in the following equations we write the charac- 
teristics Poo, Pq, V(vn) of the overlapping sphere prob- 
lem as functions of the physical parameters T/T c and 
a 3 B n c instead of r 3 n. The relation between these param- 



eters is given by Eq. (|29|). The temperature dependence 
of the spontaneous magnetization given by Eq. ( |3l| ) is 
plotted on Fig. |l2| for two experimentally viable values 
of a 3 B n c . As already mentioned before, such strongly con- 
cave M(T) behavior is often obscpfed in insulating DMS 
systems, see for example Rcfs.claEa 

The magnetic susceptibility of a sample has contribu- 
tions coming from polaron clusters and free spins. Us- 
ing the classical expression for the susceptibility of a free 



FIG. 12: Spontaneous magnetization as a function of temper- 
ature. 
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FIG. 13: Temperature dependence of the magnetic suscepti- 
bility for n c /rii = 0.02. 



magnetic moment mo 

and the fact that the total spin of a cluster with volume 
v equals vrii we arrive at 



x = h 9S 



P(tin c ;-,a B ?i c | (vn\) dv 



T 



(32) 



The first term in brackets diverges as (T — T c ) 7 with 
7 sw 1.7 at T — * T c . The second term in brackets does 
not have any singularity at T — > T c . For the tem per ature 
dependence of the susceptibility given by Eq. (|32|), see 
Fig. |l|. 

Using the classical magnetization relation 



(33) 



where C(x) = Boo{x) — cotan x — 1/x is the Langevin 
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FIG. 14: Temperature dependence of magnetization in finite 
magnetic field with Ez = g\^BSB. 
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FIG. 15: Magnetic field dependence of the sample magneti- 
zation at given temperature with Ez = g\^,BSB. 



function, wc obtain the expression for the magnetic mo- 
ment per unit volume in finite magnetic field 



M(B) 



T 



V I vn c ; — ,a%n c ] vn\ C ( ^ 1 ^' B _ 1 ] <!,< 



T 



riiSPo 



T 
T, 



,a B n c C 



gHJ, B SB \ 
T ) 



(34) 



Figures |4| and [l^ illustrate the dependence of the mag- 
netization on the temperature and magnetic field respec- 
tively. Since the spins of polarons and their clusters are 
much large than those of three impurities, the magne- 
tization curve of Fig. |l5| has two characteristic scales. 
First, polarons and their clusters are polarized (inset of 
Fig. |l5| ), then the free spins are polarized at much larger 
fields. 

We mention that our percolation-theory critical expo- 
nents are 7 w 1.7 and (3 ~ 0.4 for the susceptibility 
(Fig. |l3|) and magnetization (Fig. [l2|), as compared with 
the mean-field results of 7 = 1 and (3 — 0.5 and the best 
existing numerical estimates of 7 w 1.39 avid (3 « 0.37 
for the three-dimensional Heisenberg modelJ-^3 



III. CONCLUSION 

The temperature-dependent magnetization results pre- 
sented in this work for different but interrelated theo- 
retical models are qualitatively consistent with one an- 
other. In particular, the two theories for metallic DMS 
systems with itinerant carriers, namely the degenerate- 
carrier Weiss static mean- field theory (see Sec. [IB 2) and 
the dynamical mean-field theory (see Sec. [IC), both give 
onwardly convex M(T) carriers, similar (but not identi- 
cal) to that in the textbook molecular mean- field result, 
with the convexity being enhanced (suppressed) with in- 
creasing (decreasing) delocalized carrier density. In typi- 
cal situations involving itinerant carriers in metallic DMS 



systems, we find temperature dependent magnetization 
curves which are almost linear at lower temperatures in 
excellent agreement with experimental observations in 
metallic GaMnAs systems. Such apparent "non-mean- 
field-like" magnetization behavior, found here both in our 
delocalized-carrier static and dynamic mean-field theo- 
ries, arises from a combination of reasons: (1) the double 
Brillouin function [i. e. B(B(x))] form of the coupled field 
magnetization [see Eq. (g)]; (2) the low values of n c /n\, 
leading to the Mn moments feeling on the average a much 
reduced number of free carriers. As the carrier density 
is increased, for example by suitable annealing in recent 
experimental studies, such linear magnetization curves 
evolve toward the more conventional outwardly convex 
magnetization, as can be seen ip.puii-results and in re- 
cent annealing experiments.^E§liZlll3 : tj In the case of lo- 
calized carriers appropriate for insulating DMS systems 
(many of which are also found to be ferromagnetic with 
well-defined Curie temperatures, e. g. Qai-ajMn^As for 
x < 0.03, Ini-^Mn^ASjQ and Gei-^MnJj), our two com- 
plementary theories, namely the n on-degenerate-carrier 
static mean- field theory (Sec. II B 1 ) and the percolation 
theory (Sec. II D), give qualitatively similar magnetiza- 
tion behavior. In particular, at low carrier densities the 
M(T) curves in strongly insulating DMS systems exhibit 
strikingly "non-mean-field-like" outwardly concave mag- 
netization behavior, as often observed in insulating DMS 
systems. Again, this unusual concave magnetization be- 
havior arises from a combination of the strongly local- 
ized nature of the carrier system and the low values of 
the carrier density (i. e. n c jn\ <C 1) in DMS materials. 
We emphasize, however, a very important feature of the- 
oretical results presented in this work: the convex mag- 
netization behavior in the strongly metallic case and the 
concave behavior in the insulating case are not a sharp 
dichotomy in DMS properties — these are really the two 
extremes of a continuum of possible magnetization be- 
havior in D MS m aterials, as can be inferred from the re- 
sults of Sec. [IB. The typical DMS magnetization behav- 



ior should lie somewhere in between these two extremes of 
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highly concave (low carrier density and strongly localized 
insulating DMS) and highly convex (high carrier density 
and strongly delocalized metallic DMS) magnetization, 
leading to the generic experimental observation of almost 
linear M(T) behavior in ferromagnetic Gai-^Mn^As. 

Another often- mentioned peculiar aspect of DMS mag- 
netization, namely the low value of the saturation mag- 
netization (compared with the number of Mn ions in the 
sample as given by the value of a; in Gai-^Mn^As), is also 
apparent in our theoretical results presented in this work. 
In particular, the presence of direct antifcrromagnetic 
coupling between the Mn moments may drastically sup- 
press the saturation magnetization, as can be seen in the 
results presented in Figs. || and || of this paper. This ab- 
sence of complete magnetization saturation in our mean- 
filed results is akin to having effective ferrimagnetism in 
the system. In addition, the calculated magnetization 
is may be much lower than the saturation value except 
at very low temperatures, see Fig. E2[ which could also 
explain the lack of magnetization saturation. It is likely 
that in the real DMS materials an appreciable fraction of 
the Mn moments are magnetically inactive because they 
do not sit in cation substitutional sites, but rather at 
defect sites such as interstitials and are antiferromagnet- 
ically coupled to other Mn moments. Such magnetically 
inactive Mn atoms could be one of the reasons for the low 
values of saturation magnetization in Gai-^Mn^As. We 
also note that our percolation theory provides another 
possible explanation for the observed low values of sat- 
uration magnetization in low carrier density insulating 
DMS materials. Since the infinite cluster of percolating 
bound magnetic polarons triggering the long-range ferro- 
magnetic ordering necessarily leaves out a large number 
of Mn moments (which are not parts of the infinite clus- 
ter except at T — > 0), one naturally expects a very low 
saturation magnetization except perhaps at T -C T c . 

A recent series of potentially important annealing 
experiments in metallic Gai-^Mn^As samples by sev- 
eral different groups may eventually shed considerable 
light on our understanding of the mechanisms underlying 



DMS ferromagnetism. These experiments, while differ- 
ing somewhat on the details, all find that suitable low 
temperature "optimal" annealing may enhance the mag- 
netic properties of Gai-^Mn^As by increasing T c and 
more importantly for our purpose, by enhancing M(T) 
to more convex almost standard mean-field-like behavior. 
This enhancement of magnetization seems to correlate 
well with improvement in the metallicity of the annealed 
sample (e. g. higher conductivity) and with an increase 
in the hole density. Annealing may also be enhancing 
the magnetic properties by annealing away some of the 
Mn interstitials. Our theoretical results are in excellent 
agreement with these annealing experiments, since in- 
creasing n c jri\ does lead to enhanced (and more convex) 
magnetization in our theory. In this context, it will be 
very helpful to have more detailed information on the 
M(T) behavior in strongly localized insulating samples 
both for Gai-^Mn^As with x < 0.03 and for other insu- 
lating DMS materials (e. g. Gei_ x Mn x ). 

Our final comment addresses the role of disorder (es- 
sentially neglected in our work, except for the percolation 
theory part, as we treat it in the virtual crystal approx- 
imation, i. e., we assume that holes and local moments 
only feel the average effective field ), which is invariably 
strong (and not yet well understood) in DMS materials. 
In particular, even the "metallic" DMS systems are in 
effect very poor metals with mean-free paths which are 
at or below the Ioffe-Regel limit (with low-temperature 
mobilities of the order of a few cm 2 /Vi*iOnly), indicating 
the presence of very strong disorderly In the presence 
of such strong disorder, various spin glass ground states 
may compete with ferromagnetic ground states. We be- 
lieve that it is important to look for signatures of spin 
glass physics in low-temperature DMS magnetic prop- 
erties. It may very well turn out that spin glass phases 
dominate the regime of parameter space (e. g. x < 1% or 
x > 10% in Gai-^Mn^As) where a ferromagnetic ground 
state does not seem to stabilize in DMS materials. 

This work is supported by US-ONR and DARPA. 
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